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Abstract 

Density and temperature conditions in the solar core suggest that the micro- 
scopic diffusion of electrons and ions could be nonstandard: Diffusion and fric- 
tion coefficients are energy dependent, collisions are not two-body processes 
and retain memory beyond the single scattering event. A direct consequence 
of nonstandard diffusion is that the equilibrium energy distribution of parti- 
cles departs from the Maxwellian one (tails goes to zero more slowly or faster 
than exponentially) modifying the reaction rates. This effect is qualitatively 
different from temperature and/or composition modification: Small changes 
in the number of particles in the distribution tails can strongly modify the 
rates without affecting bulk properties, such as the sound speed or hydro- 
static equilibrium, which depend on the mean values from the distribution. 
This mechanism can considerably increase the range of predictions for the 
neutrino fluxes allowed by the current experimental values (cross sections and 
solar properties) and can be used to reduce the discrepancy between these 
predictions and the solar neutrino experiments. 
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I. INTRODUCTION 



Nowadays solar modeling seems to have reached a satisfactory stage [|T]-|^ . The inclusion 
in the latest models of higher-order effects, such as the diffusion of heavy elements, brings the 
theoretical predictions in good agreement even with the detailed helioseismological data [^|| . 

Nevertheless, solar neutrino experiments still contradict the model predictions. The 
problem appears to be almost independent of the details of solar models 0"15|, and it has 
mostly been interpreted as a hint of new physics PJT6|- p7| , p^ . However, not everybody 



believes that the discrepancy is large enough {3a effects have often disappeared) and that 
the solar models, in spite of their successes, are still solid enough to pose a real problem pS 



In this context, there has been a considerable amount of work devoted to answering 
questions such as: how large are the uncertainties of the solar model input parameters? Has 
something been left out of standard solar models? How does this affects predictions for the 
fluxes and the status of solar neutrino problem (SNP) |p5|-53. 



The solar core is a dense, strongly interacting many-body system: no realistic microscopic 
calculation of such a system exists. However, only three properties are important for the 
purpose of solar modeling. The first one is the equation of state, i.e., the relation between 
the average one-body local properties of the system (temperature, pressure, density). It is 
the equation of state that regulates the local hydrodynamical equilibrium and, in addition, 
contributes to the interpretation of the helioseismological data. The second one is the local 
opacity, which controls the energy transmission. In this paper we are interested in the last 
feature, the two-body relative-energy distribution or two-body correlation between particles. 
The rates of the most important reactions in the Sun are strongly affected by the high-energy 
tail of this distribution. 

Solar models implicitly assume that the solar core can be described in terms of a gas of 
particles interacting via two-body short-range forces with no many-body effects apart for 
mean-field screening. Therefore, the velocity distribution and, in particular, the relative 
velocity distribution for the particles involved in the reactions are Maxwellian. 



Two decades ago, Kocharov and collaborators [p4|-|36||, and Clayton and collabora- 
tors P7| , |55| speculated that the high-energy tail of the relative energy distribution could 
depart from the Maxwell-Boltzmann (MB) exponential form, exp{—E/kT}. At that time 
the SNP consisted only in the low yield of the Chlorine experiment compared to theo- 
retical predictions |^ . Since only the highest-energy proton can significantly penetrate the 
Coulomb barrier of the reaction ^Be(p, 7)'^B, Clayton suggested that a significantly lower 
boron-neutrino flux could be obtained by depleting the small number of protons in the tail 
and proposed to parameterize the small deviation with a Gaussian factor, exp{—6{E/kT)'^}: 
a value of 5 = 0.01 was sufficient to solve the SNP. This suggestion met some criticism and 
has been mostly ignored (see, however, a few recent papers |^I|-^). In particular, it has 
been remarked that deviations from the MB energy distribution could not be due to nuclear 
reactions keeping the distribution out of equilibrium |^5|-|47|fll. In addition, nobody was 
able to give a microscopical dynamical derivation of a non-Maxwellian distribution for the 
specific conditions inside the Sun (however, nobody has demonstrated with realistic calcula- 
tions, which do not already assume the validity of the Boltzmann transport equation, that 
the distribution in the Sun is Maxwellian either). 
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Several new developments have convinced us that it is necessary to reconsider the possi- 
bility of deviation from a MB distribution. First of all, non-Maxwellian equilibrium energy 
distributions have been shown theoretically possible P8|-|52| and relevant to many physical 
systems [^3H63|1. Then, it has been recently proved by explicitly solving the Fokker-Planck 
equation that a velocity-dependent diffusion or friction coefficient results in an equilibrium 



energy distribution that departs from the MB one [54,65|. One should consider that con- 



stant diffusion and friction coefficients are only the first terms of a derivative expansion of 
the linearized hydrodynamic equations: Few percent contributions from the next terms in 
the expansion are not unrealistic, and are sufficient to give appreciable deviations from the 
MB distribution. Moreover, successive scattering events should not be independent (Marko- 
vian) in the solar interior, but correlated over periods corresponding to a few scattering 
processes. Indeed, such correlations have also been found to be related to anomalous dif- 
fusion [^,^ and to the Tsallis' statistics that is a generalization of the Boltzmann-Gibbs 
statistics Time correlations, velocity dependence of the transport coefficients or 

other deviations from the so-called Navier-Stokes limit of the hydrodynamic equations are 
all manifestations of multiple particle collisions and of the many-body nature of the system, 
which is not completely described by two-body short-range interactions. In addition, the new 
data of the solar neutrino experiments (Chlorine [^], Kamiokande pU| , [7n[ , Gallex [[7I|-[73|] 
and SAGE |]73|) have made the actual SNP more puzzling |5,75] and the possibility that 
Clayton's suggestion could contribute to its solution should be reassessed. 

It should be stressed from the beginning that modifications of the shape of the energy 
distribution are not equivalent for solar models to changes of temperature and/or densities; 
rather, a new degree of freedom is introduced. For instance, one might think that the effects 
of a distribution with a depleted high-energy tail could be reproduced by a MB distribution 
with a lower temperature. However, a lower temperature produces two effects: (1) the 
rates are reduced by changing the thermal average (va) (most of the contribution to this 
average comes from the high-energy tail because of the particular energy dependence of 
the main cross sections a); (2) the system finds a new hydrostatic equilibrium because the 
average momentum of the particles becomes smaller (particles of all energies give comparable 
contributions to this average). In contrast, a change of shape of the distribution can reduce 
the rates and, at the same time, maintain the same hydrostatic equilibrium, since the two 
effects are dominated by particles from different parts of the energy spectrum. Indeed, the 
effects on solar models of changing the shape of the energy distribution could be reproduced 
by simultaneous local changes of temperature and cross sections. 

Given the widespread misconception about the inevitability of the Maxwell-Boltzmann 
velocity distribution in the solar core, we briefly review in Sect, y the general assump- 
tions under which the MB distribution is derived, show that these assumptions can be only 
approximately verified in the solar interior, and discuss some of the expected corrections 
to the standard treatment. In particular we shall consider (Sect. |IIID three concrete ex- 
amples of anomalous diffusion (velocity-dependent corrections to the diffusion coefficient, to 
the friction coefficient, and slowly-decaying velocity autocorrelation), and how they generate 
deviations from the MB energy distribution. Then in Sect. I^we illustrate the consequences 
of small deviations for solar neutrino fluxes and estimate the magnitude of the deviations 
necessary to change the fluxes by amounts relevant for the SNP. We reserve Sect. to our 
conclusions. 
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II. THE MAXWELL-BOLTZMANN DISTRIBUTION AND THE SOLAR CORE 



In this Section we briefly review the standard hypotheses that lead to the MB distribution 
for the single-particle velocity, and the ones under which the two-particle relative-velocity 
distribution is also Maxwellian. Then we estimate the order of magnitude of the relevant 
physical quantities and show that these hypotheses are not met, making it likely that the 
MB distribution be only an approximation to the real distribution. These same "order of 
magnitude" considerations suggest some of the corrections to the hydrodynamical equations. 

A. The ubiquity of the Maxwell-Boltzmann distribution 

The velocity distribution can be studied either with a dynamical kinetical approach or 
with the methods of equilibrium statistical mechanics [|76|-|80[|. 

The kinetical approach yields a hierarchy of differential equations (the Bogoliubov-Born- 
Green-Kirkwood-Yvon hierarchy). The n-th differential equation is an equation of mo- 
tion for the n-body distribution function /„(ri, f i; . . . ; r„, f„; t), consisting of a streaming 
term, which involves /„, and a collision term, which involves the (n + l)-body distribu- 
tion function. In particular, the one-body distribution function /i(ri,t>i;t) verifies a dif- 
ferential equation that involves the two-body distribution function /2('"i, "^i; r2, f2; t) in the 
collision integral. The derivation of the BBGKY hierarchy already assumes that (1) the 
collision time is much smaller than the mean time between collisions and that (2) the in- 
teraction is sufficiently local. When these conditions are not met, collisions are not well- 
defined events in space and time. The rate of change of the /'s at time t would depend 
not only on the /'s themselves at time t, but also on its previous history (the process 
is not Markovian), and higher-order spatial derivative cannot be discharged. Assuming 
that the BBGKY hierarchy is valid, the additional hypothesis of molecular chaos (Boltz- 
mann's Stosszahlansatz), i.e. that (3) the velocities of two particles at the same point are 
not correlated /2(r, t>i; r, t>2; t) = fi{r,vi]t) x /i(r, t>2;t), allows the truncation of the hier- 
archy of equations and yields the Boltzmann transport equation for the one-body distri- 
bution function fi{r,V]t) alone. From the Boltzmann transport equation the Boltzmann's 
H theorem can be derived whose consequences are that (i) under arbitrary initial condi- 
tions limf^oo f{v,t) = /o(f) Q, and that (ii) /o is the equilibrium distribution if and only if 
fo{vi)fo{v2) = fo{v[)fo{v2), where (fi,f2) and {v[,V2) are the velocities before and after the 
collision, i.e. J2i^og fo{vi) is conserved in the collision. Then the additional assumption that 
(4) energy is locally conserved when using only the degrees of freedom of the colliding par- 
ticles yields that the energy of the particle contributes linearly to log/o(f), and, therefore, 
the one-particle equilibrium distribution is Maxwellian. 

The equilibrium statistical mechanics approach uses the concept of most probable value 
of the distribution, and the large number of degrees of freedom assures that large fluctuations 
away from the most probable distribution are extremely unlikely in non-critical conditions. 



For the sake of this discussion, we drop from now on the spatial label r, that should be reintro- 
duced when discussing non-local interactions or density gradients 
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Probabilities are calculated using the Postulate of Equal a Priori Probability and, when 
deriving the MB distribution, the assumption (3) that the velocity probabilities of different 
particles are independent. In addition, it is necessary that the total energy of the system 
could be expressed as a sum of a term quadratic in the momentum of the particle and 
independent of the other variables, and a term independent of momentum (this second 
term includes the energy of the rest of the system and the interaction energy). This second 
condition is related not only to assumption but also to (1) and (2), since (1) and (2) 
are not verified the resulting effective two-body interaction is not local and depends on the 
momentum and energy of the particles. 

Finally, even when the one-particle distribution is Maxwellian, assumptions (3) and (4), 
i.e. that the velocities are uncorrelated and energy is locally conserved, are again necessary 
to deduce that the relative-velocity distribution is also Maxwellian. 

At last, we suggest two of the reasons of the ubiquity of the MB distribution. On the one 
side, whenever the above assumptions are a good approximation (typical examples are those 
systems that are dilute in the appropriate variables), the resulting equilibrium distribution 
can be demonstrated to be Maxwellian independently of the details of the interaction. On 
the other side, even when conditions are such that one can not rigorously deduce the form 
of the distribution, considering that most of the experimental measurable observables do 
not really test the form of the distribution, but only a few of its moments, the assumption 
of a MB form has several advantages: (a) It becomes exact in the dilute limit; (b) Being 
determined by a single parameter or scale, the second moment of the velocity distribution, 
it is consistent with our maximal ignorance (it assumes the least about the distribution); 
(c) It allows simple analytical treatments. 

B. The solar core quasi-plasma 

In the light of the above considerations, we examine the situation in the solar core, and 
discuss what kind of corrections to the simplified picture of a dilute gas could better describe 
a system where, as we shall argue, multiple particle collisions and many-body correlations 
between clusters of tens of particles should be present. Of course, our discussion suggests 
only qualitative answers, since these issues need a more sophisticated theoretical framework, 
better microscopical models and, perhaps, they could still be quantitatively settled only by 
a numerical dynamical microscopic calculation, which is not trivial and, to our knowledge, 
still missing. 

Ions in the solar core are completely ionized: T^, = 1.36 keV. 

Only the many-body correlations between ions and electrons cut off the range of the 
electromagnetic interaction. The Debye-Hiickel estimate (weak mean-field static approach) 
gives a screening length {R\) = kT / {ATre"^ J2i ^ii^i)) of the order of the average interparticle 
distance in the center of the Sun: Rd cm, where n is the average density. 

In addition, the same Debye-Hiickel approach estimates that, at the average interparticle 
distance, i.e., when the potential energy is at its minimum, the potential energy is already 
about 4% of the kinetic energy: its contribution is comparable to the kinetic energy for large 
portions of the particle trajectories. 

These "back of the envelope" calculations already suggest that it is not possible to 
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completely describe the system as streaming particles plus collisions. Since the collision time 
is of the same order of the mean time between collisions and the size of the quasi-particle 
(ion plus screening cloud) is of the order of the distance between particles, a description in 
terms of a gas of compact quasi-particles interacting via local two-body collisions can only 
be an approximation whose corrections involve higher order derivative terms. In fact, the 
relevant expansion parameter in the BBGKY hierarchy is nrj], where n is the density and 
ro the range of the interaction; but uTq ^ 1 using ro ~ Rd- 

Moreover, while the Debye-Hiickel screening strictly applies only to a charge at rest 
or travelling at constant velocity, i.e. in the limit that the plasma reaction time is in- 
finitely fast compared to the particle rate of velocity change, the real screening is dy- 
namical and the inverse plasma frequency gives us an estimate of the reaction time: 
tpi = Lj~i^ = \Jm/ (47rne^) ~ 10"^''' sec. Times of this order of magnitude are necessary 
to build up screening after a hard collision and should be compared to the typical collision 
time tcou = {(Jvn)~^ ~ 10~^^ sec. 

We describe the fact that tpi ^ tcou by saying that collective effects have a time scale 
comparable to the average time between collisions and, therefore, several collisions are neces- 
sary before the particle looses memory of the initial state. Another consequence of tpi ~ tcoii 
is that collisions between the quasi-particles (bare particle plus the screening cloud) are in- 
elastic: part of the energy is dissipated by the process of the creation and successive removal 
of the screening cloud. 

Since the plasma parameter, i.e., the number of particles in the Debye sphere, is not 
large {Njj = {4:7i/3)R^n ^ 1), a description in terms of a high-temperature plasma {1/Nd ~ 
7^(1/2) 2^- (3/2) _^ jg j^Q^ ^ good approximation either. We could name it a quasi-plasma: 
strong many-body effects on the scale of a few average interparticle distances appear to be 
necessary, and clusters of a few (or tens) of correlated particle participate in the collisions, 
the number depending on the time scale considered. 

In the next Section, we show that when the lowest-order kinetic equations, which are 
valid in the low-density limit, are supplemented with terms that have the form suggested 
by the above considerations, i.e. higher derivative terms (spatial nonlocality), time corre- 
lations (temporal nonlocality), or nonlinearity, the resulting equilibrium distribution is not 
Maxwellian. 



III. ANOMALOUS DIFFUSION, TIME CORRELATIONS AND NONSTANDARD 

STATISTICS 

For concreteness we consider three possible corrections: i) A correction to the lowest order 
friction coefficient J(f ); ii) a correction to the lowest order diffusion coefficient D{v)\ and iii) 
a modification of the two-body time correlations. These three types of corrections affect the 
Maxwell-Boltzmann statistical distribution. Indeed, similar corrections have already been 
shown to exist in hydrodynamic systems [^,^,0: direct microscopic calculations could 



prove this possibility also in the solar context. 

It is clear that such corrections imply other consequences. For instance, it is well-known 
that the present approach to the slow diffusion of heavy elements could also be modi- 
fied In the present paper, we are interested in only one of these consequences: the 
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actual equilibrium statistical distribution of the relative-energy departs from the Maxwell- 
Boltzmann equilibrium distribution. 

Corrections i) and ii) have already been considered in a more general context p^ , |65| , 
and we only recall the main points. We assume that the system is not too far from the 
standard regime that leads to the MB distribution, so that an expansion starting from the 
usual formalism makes sense. The Fokker-Planck equation, given in the Landau form, is 



d_ 
di 



f{t,v) 



d_ 

dv 



J{v)f{t,v) 



d_ 

dv 



D{v)f{t,v) 



(1) 



where /(t, v) is the distribution probability of particles with velocity v at time t and J{v) and 
D{v) are the dynamical friction and diffusion coefficients. The stationary distributions are 
the asjTiiptotic solutions of the above equation. To lowest order J{y) = v/t and D{v) = e/r, 
where the constant r > has dimension of time (m/r is the friction constant) and y/e has 
dimension of a velocity (e = kT/m for Brownian motion). At equilibrium one obtains the 
well-known Maxwellian distribution 



f{v) = lim f{t,v) 

t — >OQ 



exp 



mv 
'2¥f 



(2) 



We can generalize the standard Brownian kinetics considering the expressions of the 
quantities J{v) and D{y) to the next order in the velocity variable: J{v) = v/t {1 + Piv'^/e) 
and -D(f) = e/r (1 + 7it>^/e); these higher derivative terms can be interpreted as signals of 
nonlocality in the Fokker-Planck equation. 

If /5i = and 71 7^ we find the Tsallis' distribution 



fiv) 



1 + (g - 1 



mv 



1/(1-9) 



e 



1 + (g - 1) 



mv 
2kf 



(3) 



where g — 1 = 271/(271 -|- 1), is the Heaviside step-function, and kT/m = e(2 — q). When 
the characteristic parameter q is smaller that 1 (—1/2 < 71 < 0), this distribution has a 
upper cut-off: mv'^/2 < kT / (1 — q) (the tail is depleted). The distribution correctly reduces 
to the exponential Maxwell-Boltzmann distribution in the limit g — > 1 (71 0). When 
the parameter q is greater than 1 (71 > 0), there is no cut-off and the (power-law) decay is 
slower than exponential (the tail is enhanced). 

If /3i 7^ and 71 = 0, we find a Druyvenstein-like distribution: 



/(-)~exp<l-^-/3i 



(4) 



which has also the functional form suggested by Clayton to parameterize a small deviation 
(depletion) from the Maxwellian statistics. 

The statistical distribution of Eq. (0) has an additional appealing feature: it naturally 
appears in the context of the generalized Boltzmann-Gibbs statistics obtained by introducing 
a new non-extensive entropy (Tsallis' entropy): 



k 



('?-l)^ 



(5) 
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The formal structure of the conventional thermostatistics is maintained and its results are 
also naturally generalized 1 48, 56. 5^. Apart from the formal aspect, this distribution is also 
attractive because of the many systems where it plays a role []53|-|63[1. Non-extensivity (g 7^ 1) 
arises in systems with long-range interactions (gravitational systems []63 
condensed matter 



plasmas ||85| ,|^ 



B/]) or with long memory at the microscopic level. 



We mention also the possible connection between time correlations and extended statis- 
tics (other mechanisms are also possible). For a Markovian scattering process, the time 
correlation between particle velocities is by definition proportional to a delta function in 
time: {v{0)v{t)) ~ 6{t). As discussed above, we expect that particles loose memory of the 
initial state only after a few scattering processes; we can model the long-time asymptotical 
behavior of the velocity-correlation as {v{0)v(t)) ~ t^^^^'^\ If 7 > 1, i.e., the correlation de- 
cays sufficiently fast, the diffusion process is no qualitative different from the delta-function 
case: {x'^{t)) ~ t. The same standard result holds if < 7 < 1 and < J {v{0)v{t)) < 00. 
However, it has been shown that, if < 7 < 1 and / {v{0)v{t)) = (or very small), or 

if — 1 < 7 < 0, the diffusion is anomalous (x^(t)) ~ t^+'^ (~ tlogt, if 7 = 0). Indeed, Tsal- 
lis ||5^ shows that the generalized entropy Sq quite naturally generate anomalous diffusion 
{{x'^{t)) ~ t^^"^): this same generalized entropy leads also to the non-Maxwellian probability 
distribution for the velocities given by Eq. (RT). 



IV. NONSTANDARD STATISTICS AND SOLAR NEUTRINOS 

From the considerations above we infer that generalized distributions, among which 
Tsallis' distribution has a special theoretical appeal, could better approximate the situation 
in the solar interior. 

One could proceed and study the effect of using generalized distributions on solar models, 
distributions with both depleted and enhanced tails, for small and large deviations from the 
Maxwellian statistics. In particular, it would be interesting to study how the relative balance 
among the different reaction chains would change because of non-Maxwellian statistics; 
for instance, a distribution with an enhanced high-energy tail could make the CNO cycle 
important at relatively lower temperatures. This study would be very interesting, since it 
could experimentally constrain deviations from the MB statistics, given the high sensibility 
of these reaction rates to the tail of the distribution. However, such a systematic study, 
which should be performed by consistently including nonstandard statistics in solar model 
calculations, is not the purpose of this paper. 

We shall only consider small deviations from the Maxwellian distribution and, for the 
purpose of illustration, use Clayton's parameterization with the factor Q-^i^/'^'^) . For in- 
stance, the Tsallis' distribution can be also approximated to first order in (l — q) by Clayton's 
form with 6 = {1 — q)/2 and a renormalized temperature T' = T + T(l — g). The usual 
asymptotic expansion of the integrand over the velocity distribution around the most 
effective energy, Eq, yields an analytical expression for the rate changes. This analytical ex- 
pression is valid also for 5 < as an asymptotic expansion around the 5 = case, in spite of 
the fact that the distribution is unbounded at high energy. Instead, a numerical integration 
should be performed with a suitable cut-off. In fact, the integrand decays exponentially after 
the Gamow peak when 6 > 0, and a sufficiently large cut-off does not changes the numerical 
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value of the integral. When S > 0, the integral still decays after the Gamow peak, but it 
goes back up at energies E ^ Eq/6: this contribution to the enhanced tail is "unphysical" 
and comes from the choice of the parameterization. If 5 -C 1, the large window between Eq 
and Eq/6 allows the unambiguous elimination of this contribution. The alternative is to use 
a distribution that has an enhanced tail but still decays at high energy, such as the Tsallis' 
distribution for g > 1. For small deviations from the MB distribution, the two descriptions 
give the same numerical results with the appropriate reparameterization, since one single 
number characterizes the deviation to first order. 

If one computes the thermal average (va) with the modified distribution for a two-body 
reaction with Coulomb barrier, one finds that to the leading order in 6 

= g-<5(S«/feT)2 ^ ^ (6) 

{vai)o 

where Eq is the most effective energy (maximum at the Gamow peak) ^ 

which depends on the reaction i through the charges Z and weights A of the ions, and on the 
relevant average temperature T. Here, Tc = 1.36 keV is the temperature at the center of the 
Sun. In Table | we report the values of 7^ = {Eq ^ /kT)'^ for the five most relevant reactions 
in the Sun: p + p{i = 1.1), p+^Be (z = 1.7), p+^^N {i = 1.14), ^He + ^He {i = 3.3) and ^He 
+ ^He {i = 3.4). Changing (va) for the zth reaction will affect the whole solar model and, 
in general, all fluxes will change. We estimate the effect on the fluxes by using power-law 
dependences 

R.^l^ = rj( MA = e- '^^^-^^ (8) 



for the fluxes j = ''Be, ^B, ^^N and ^^O, while we have used the solar luminosity con- 
straint [|r§] to determine the pp flux, Rpp = 1 + 0.087 x (1 - R^c) + 0.010 x (1 - R^) + 
0.009 X (1 — Rq), and kept fixed the ratio ^ = ^pep/^pp = 2.36 x 10~^. The exponents 
aij = dln^j/d ln{vai) (see Table |Tp have been taken from Ref. []T3| , where it is also discussed 



why solar models depend on (f (7)33 and (f (7)34 only through the combination (f (7)34/ ^J (f (7)33 
and why it is a good approximation to keep the ratio ^ constant. 

In principle 6 should be determined by a direct calculation of the complex many-body 
system and could be different for every reaction {6 — > 6i). The energy distribution can 
be influenced by the specific properties of the ion (charge and mass) and by the different 
conditions of the environment in those parts of the Sun where each of the reactions mostly 
takes place. However, a direct calculation is not simple and it does not exist for the solar 
interior. Therefore, for the only purpose of estimating the potential effect of nonstandard 
distributions, we consider two simple models and use the corresponding 5('s) as free param- 
eter (s). The first model assumes the same deviation S for all distributions, while the second 
model assumes that only the p + ^Be relative-energy distribution (parameterized by 6Be) and 
the two helium reactions (parameterized by Sne) are non standard. 
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In the first case, one finds by substituting Eq. (Q) in Eq. that 



^(0) ^ ' 

3 

where jSj = J2i c^jili are reported in Table |I|. This dependence of the fiuxes on S is in good 
agreement with Clayton's numerical calculation Using the model of Ref. ^ as reference 
model and the latest experimental results (see Table pT]) , we obtain the best fit for 6 = 0.005 
with a = 35. 

In the second case, we proceed similarly, but we use ^Be for the reaction p+^Be and S^e for 
the two reactions He + He: the corresponding = aji'ji\i=ij and = «j^3. 4(73. 4 — 73.3/2) 
are also reported in Table 0. As shown in Table [TTTI the best fit is obtained for ^Be = —0.018 
(negative 6 corresponds to an enhanced tail, g > 1 in Tsallis' distribution) and d^e = 0.030 
with a = 20. 

From previous analyses we already knew that it is not possible to obtain a very good fit 
to all experiments even when the fiuxes are used as free parameters |pT| , p5| - pil|JT3[ ] : a fit that 
has a low probability could only be obtained by reducing the boron fiux by a factor of two 
and the beryllium and CNO fiuxes as much as possible. Therefore, we are not surprised that 
we have not been able to obtain good fits, however we have been able to give a specific and 
physically motivated mechanism that greatly reduces the discrepancy between theory and 
experiment (the SSM has a = 74). Another physical mechanism that produces similar 
results is by introducing arbitrary screening factors |jl5|. It is also consistent that the second 
case, which produces the best fit, results in a depletion of the energy tail of the ^He and 
^He ions, so that ^Be and '^B fiuxes are strongly suppressed, and that the energy tail of the 
p+^Be reaction is enhanced so to bring up the ^B fiux towards the measured value. 

We do not claim that this result is a solution to the SNP, in the sense of providing 
a model to fit the experimental results within one (a few) sigma. Our point is only that 
deviations from standard statistics corresponding to values of S of about 1% can change the 
neutrino fiuxes of factors comparable to those that constitute the SNP. Such values of S, or 
even larger values, cannot be excluded by the present knowledge of the strong interacting 
quasi-plasma in the solar interior. It could turn out that the actual values of the neutrino 
fiuxes coming out of the Sun could result from the interplay of several mechanisms that are 



disregarded in the standard picture |28 



In the light of the above considerations, the uncertainties of the neutrino fiuxes are con- 
siderably underestimated by not considering the possibility of non-extensive distributions. 

Finally, we wish to comment on the fact that limits on the reaction rates that come 
from determinations of the sound speed through helioseismological measurements do not 
automatically apply to changes of the rates through the present mechanism. In fact, given 
the cross sections, the reaction rates change because the densities and/or the thermal aver- 
ages, (cTf), change. However, if the statistics is not changed, the thermal averages change 
when the temperature changes, and changes of temperature and/or density clearly affect the 
structure of the solar model and the sound speed. In contrast, nonstandard energy distri- 
butions make it possible to change the thermal averages (at least for those reactions whose 
main contribution comes from the high-energy tail) without affecting the properties that de- 
pend on the bulk of the distribution, such as the sound speed and/or the equation of state. 
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Therefore, non standard velocity distributions affect the hehoseismological measurements 
only insomuch as the consequent changes of the rates modify the solar structure. 

V. CONCLUSIONS 

New developments in generalized statistics combined with several qualitative hints that 
the standard approach to the solar interior is only a first approximation to the real situation, 
make it worthwhile to reconsider the early suggestion by Clayton that the energy distribution 
in the Sun could depart from the Maxwell distribution. 

In particular, wc recall that: 

(1) The conditions in the solar core (density and temperature) do not satisfy those 
requirements that would guarantee standard diffusion and Maxwell-Boltzmann velocity dis- 
tribution. 

(2) Nonstandard diffusion is most likely present: Lowest order dynamical friction and dif- 
fusion coefficients are not sufficient: the diffusion mechanism is not described by a Markovian 
chain of independent two-body scattering events and correlations persist for time intervals 
longer than the mean time of one scattering process (memory effect). 

(3) The direct consequence of these corrections to standard diffusion is that the equilib- 
rium energy distributions of electrons and ions are not Maxwellian. In particular, the tails 
of the distributions are not exponential and, therefore, the small number of particles that 
have energies large enough to participate in those reactions that are hindered by Coulomb 
barrier can be much less (more) than the one expected in the standard distribution. 

(4) Tsallis' statistics should also be considered. 

(5) Non-Maxwellian energy distributions modify the reaction rates. This phenomenon 
has the potential of increasing the range of possible values of the reaction rates well beyond 
the ones allowed by the uncertainties in the corresponding cross sections. 

(6) Unlike a change in the temperature, which has a direct effect on the hydrostatic 
equilibrium and on the the sound speed, modifications of the distribution that affect only 
the high-energy tail do not change the solar model and the sound speed: the range of neutrino 
fluxes from models that verify hehoseismological constraints could also be increased. 

(7) If one modifies the standard distribution by a Clayton's factor e~^^^/^'^^ with 6 of 
the order of 1% (such modification cannot be excluded by the present knowledge of the 
strong interacting quasi-plasma in the solar interior) the neutrino fluxes change of amounts 
comparable to those that constitute the solar neutrino problem, even if it is not possible to 
solve the SNP by only modifying the energy distributions. 
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TABLES 



TABLE L Most effective energies for thermonuclear reactions and exponents 7 that character- 
ize the change of the thermal average {va) to the leading order in (5, when the energy distribution 
changes by a factor ex.-p{—5{E / kT)"^}: {va)s = (w(t)o exp{— 57}. 

reaction Eo/kT = {Eo/kTf 



{va)ii: p + p^^R + e^ + i^ 4.8 23 

{va)ir. j9 + 7Be^8B + 7 13.8 190 

(w)33: 3jje + 3He ^ a + 2p 16.8 281 

(w)34: 3jje + 4He ^ _^ ^ I7 4 393 

{va)i^u: p + 14N^150 + 7 20.2 407 



TABLE IL The first four rows show Oij = d\n^j/dln{va)i, the logarithmic partial derivative 
of neutrino fluxes with respect to the parameter shown at the left of the row. These numbers 



are discussed in Ref. [15|. The last three rows show Pj, and Pj^'^, the logarithmic partial 



derivative of the fluxes with respect to the parameters (5's; as discussed in the text, they are linear 





■^Be 


8B 


CNO 


{va)ii 


-1.0 


-2.7 


-2.7 


{vcr)34./^/{va)33 


+0.86 


+0.92 


-0.04 


{va)n 





1 













1 




117 


277 


338.5 


fir 





190 





/3f 


140 


150 


-6.5 
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TABLE III. The first three columns show the predicted fluxes, and the predicted gallium and 
chlorine signals in the SSM Q and in the two models with nonstandard distribution described in 
the text. The last column shows the present experimental results. For the three models is also 
given the resulting by the comparison with the experimental data. 

Models 







SSM case I 
[5 = 0) {5 = 0.005) ((5Be = 


case II 
-0.018, fee = 0.030) 


Experiment 


[lO^cm ^s 

<I>15o 

[lO'^cm-^s- 




59.1 
5.15 
0.62 
0.55 

6.62 


62.2 
2.87 
0.11 
0.10 

1.65 


63.7 
0.08 
0.75 
0.67 

2.25 


2.55 ±0.21 


[SNU] 
gallium 
chlorine 




137.0 
9.3 


100 
2.84 


97 
3.34 


75±5 
2.54 ±0.20 " 






74 


35 


20 




^Weighted 
'^Weighted 


avera^ 
averaj 


;e of 2.80 ± 0.38 
?e of 76 ±8 |73| a 


69| and 2.44 ± 0.26 @ 
nd 72± 13 [74| 



=Ref. HI 
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